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HEAT TRANSFER WITH VARIABLE FLUID PROPERTIES FOR 
PRANDTL NUMBER OF 1 
By Robert G. Deissler 


SUMMARY 

, Equations were derived for the prediction of radial distributions 
of velocity and temperature for fully developed turbulent flow of gases 
in smooth tubes for the case where the variation of fluid properties 
due to temperature variation across the tube is considered. The 
equations apply for both heat addition to the gas and heat extraction 
or. cooling.' The effect of frictional heating on temperature distri- 
bution was investigated and the analysis indicated that the velocity 
and total -temperature profiles are similar for a fluid with a Prandtl 
number of 1. 

By use of the equations for velocity and temperature distributions, 
relations were obtained among Nusselt number, Reynolds number, and 
friction factor for a fluid with a Prandtl .number of 1. It was fo und 
that the effect of variation of gas properties across the tube on the 
Nusselt number and friction-factor correlations can be eliminated by 
evaluating the gas properties at a temperature' close to the average of 
the wall and bulk temperatures. The Nusselt numbers and friction factors 
analytically predicted agree with those experimentally determined by 
measurement of average surface heat -transfer coefficients and friction 
factors in tubes for Reynolds numbers above 10,000. 


INTRODUCTION 

Although some work has been done on laminar flow with variable 
fluid properties, very little has been done on turbulent flow with 
variable fluid properties. In reference 1, equations are developed for 
velocity distributions for fully developed adiabatic turbulent flow in 
smooth tubes . The plots of the equations agree closely with experimental 
data for adiabatic turbulent velocity distributions when appropriate 
values of the experimental constants appearing in the equations are used. 
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The investigation reported herein, which was conducted at the 
NACA Lewis laboratory, is an extension of the analysis developed in 
reference 1 for adiabatic flow in smooth tubes to include flow with 
heat transfer. For flow with high rates of heat transfer, the temper- 
ature variation, and thus the variation of fluid properties across the 
tube, may be so large that it might be expected that the velocity 
distributions for flow with heat transfer would differ from those for 
flow without heat transfer. In developing the equations for turbulent 
flow of gases with heat transfer where the variation of gas properties 
is considered, the same hypotheses for the turbulent transfer of momentum 
are used as were found to apply to flow without heat transfer. Because 
of these and other assumptions made in the analysis, the equations for 
velocity distributions with heat transfer should be experimentally 
verified. 


ANALYSIS 

Turbulent Transfer of Heat 

When turbulent flow exists through a tube, portions of the fluid 
move about in random fashion (reference 1) . If heat transfer between 
the tube wall and the fluid takes place, a temperature gradient, occurs 
- across the tube and some of the portions of fluid move transversely into 
regions of different mean temperature. This motion produces heat trans- 
fer in addition to that which takes place by local molecular conduction. 

By analogy with the law for transfer of heat by molecular conduction 
q v = - k dt/dy, the equation for the transfer of heat by turbulent motion 

is often written in the following form (All symbols are defined in the • 
appendix.)* ' 

,1 t = -Pgcp6 h dt/dy (1) 

I 

where pgc^e^ is comparable to the conductivity k for molecular heat 
transfer and e h . is the coefficient of eddy diffusivity for heat with a 
value determined by the amount and the kind of turbulent mixing at a 
point. Adding the molecular heat transfer to the turbulent heat transfer 
gives for the total heat transfer^ 

q = .- (k + pgc p £ h ) dt/dy ' (2) 

As in the case of momentum transfer, the analogy between molecular 
and turbulent heat transfer is inexact because the mechanism for heat 
transfer by molecular motion is somewhat different than the mechanism 
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for heat transfer when portions of the fluid move about in random 
fashion. In the case of heat transfer by molecular conduction, transfer 
takes place only at the instant the molecules collide; whereas in the 
case of turbulent heat transfer, the transfer can take place continu- 
ously by convection inasmuch as the portions of fluid are continuously 
in contact. This difference can, however, be absorbed in the value of 
€jj, the value for which depends on the turbulence mechanism, so that 

equation (l) should still apply. 


Velocity-Distribution Equations for Flow of 
Gases with Heat Transfer 


The equation for heat transfer (equation (2)) can be rewritten as 
follows : 



The equation for shear stress (reference l) is 

T = (M PC) (4) 

Assumptions . - The following assumptions are made in the use of 
equations (3) and (4) for obtaining velocity distributions with heat 
transfer : 

1. The eddy diffusivities for momentum and heat transfer (c and 
C^) are equal. Previous analyses for flow in tubes based on this 
assumption yielded heat-transfer coefficients and friction factors that 
agree with experiment (reference 2). 

2. The expressions for eddy diffusivity that are found in reference 1 
to apply to flow without heat transfer apply also to flow with heat 
transfer with variable fluid properties. These expressions are 


€ 


n 2 uy 


( 5 ) 


for flow close to the wall and the' Karmen relation 
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for flow far from the wall, where n and k are experimental constants 
having the values 0.109 and 0.36, respectively. 

3 . The ,shear stress t and the heat transfer per unit area q 
can he considered' uniform across the tube. It is shown in reference 1 
that the shear stress can he considered, with good approximation, to be 
uniform across the tube for turbulent flow. In the present analysis, 
the heat transfer per unit surface area q is also considered uniform 
across the tube. 

: 

4. The static pressure can be considered constant across the tube. 

5. The Erandtl number (Pr = l) and specific heat can be considered 
constant with temperature variation. 


6-. The effects of frictional heating (compressibility effects) on 
velocity distributions can be neglected- for subsonic flow (reference 1) . 


Halations between gas -property values and velocity . - When heat 
transfer takes places in a ga3 flowing through a tube, a temperature 
gradient exists across the tube and the gas properties vary from point to 
point. The gas properties that vary appreciably with temperature and are 
important in the present investigation are the density p, the viscos- 
ity p, and the thermal conductivity k. The specific heat is con- 
sidered constant because its variation with temperature is small compared 
with that of the other properties. 

Dividing equation (3) by equation (4) then gives, for a Prandtl 
number of 1 and constant shear stress and heat transfer across the tube, 


*0 

c p^0 


dt 

du 


( 7 ) 


Integration of this equation gives 
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^0 

c p« T 0 


du- = 


U 


0 



J 0 


4 o . . 

— _ — u = t n - t 

c p 8 T 0 ° 


( 8 ) 


This equation can he made dimensionless by substitution of the dimension- 
less quantities 



and 



The equation then becomes 


3“ = 1 - Pu + • (9) 

t o 

where 3 is the heat -transfer parameter. Equation (9) indicates that . 
the velocity and temperature profiles are similar for a Prandtl number. 
of 1 so that if the velocity profile is known, the temperature profile 
is also known. 

The variation of density across the tube can then be written by 
using the assumption of constant static pressure across the tube and 
the perfect gas law as 


_p_ = ^0 = 1 

Pq t 1 " pu + 
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or the density at any point across the tube is 



( 10 ) 


The variation of the viscosity of a 
shown from viscosity data to he given hy 


gas with temperature can he 
an equation of the form 



The viscosity at any point across the tube is then 

P = Pq (- 1 " 3 u+ ) d 


( 11 ) 


( 12 ) 


With the assumptions of constant Prandtl number and constant specific 
heat, the law for variation with temperature or velocity for the conduc- 
tivity must he the same as that for the viscosity. 

Flow dose to wall . - For obtaining the equation for flow close to 
a smooth wall, the expression for e from equation (5) is substituted 
into equation (4) to give 


t q ■ (p + puy n 2 ) (13) 

where t is replaced by t q because the shear stress is assumed not 
to vary with radius. In order to obtain equations that are in dimension- 
less form, the following well-known dimensionless quantities are used: 


‘ V t o/po 

+ . V t <Ao 

y “ Pq/ p 0 


(14) 


y 


(15) 
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Then 


and 


du _ du + 
dy " n Q dy + 


(16) 


a 2 u _ T oV T o/ p o d 2 u + 

dy 2 t^0 2 / p 0 dy +2 


(17) ' 


Substituting the expressions for density and viscosity from equations (10)\ 
and (12) and also equations (14) to (16) into equation (13) in order to 
make it dimensionless gives 


n 2 u + y + du + 
1 - 0u + dy + 


(1 - 3u + ) 


d du^_ 
dy 1 " 


By rearrangement , 


d£ . nVii . (1 . p u+) a 

du + l-pu + 


(18) 


Equation (18) is a first-order linear differential equation, the 
solution for which is, with the condition that u + = 0 ,when y + = 


nu+ 


n 2 u + 


n 


y 


= e 6 (l - pu’ 1 ’) 


n 2 u + 


(1 - 


n 2 , 

“P + d 

pu + ) P du + 


0, 


(19) 


m 


For flow close td a wall with heat transfer, equation (19) gives the 
relation between u + and y + for various values of the heat -transfer 
parameter p . 

Flow at a distance from wall . - For obtaining the equation for 
flow at a distance from a wall, the viscous shear stress is- assumed 
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negligible compared with the turbulent shear stress. Substituting the 
expression for e from equation (6) into equation (4) and substituting 
t 0 for t (assumption 3) gives 



/ . 

Substituting equations (10), (16), and (17.) into equation (20) gives 


d 2 u + m _ k ( tatf 

dy +2 Vl - pu+ WV 


where the minus sign was chosen in taking the square root in order to 
make * positive. After one integration, equation (21) becomes 



Equation (23) can be integrated to yield 



Cjp 

2k 2 




pu' 



+ C 


As the wall is approached, the velocity gradient becomes very large 
compared with that at a distance from the wall so that dy + /du + approaches 
zero as y+ approaches zero. By substituting dy+/du + from equation (22) 
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“J- •fi ~ P u+ / 

for e in the preceding equation and letting dy + /du + 

equal zero when y + and u + equal zero, C equals zero so that 

‘ $ 6 " ♦ x ) <* 

For flow at a distance from a wall with heat transfer, equation (24) 
gives the relation between u+ and y + for various values of the heat- 
transfer parameter 3. For 3 = 0, equation (24) is indeterminate. For 
this case (3 is set equal to zero in equation (21) before integrating, 
and the well-known logarithmic equation is obtained: 

/ 

u + = ^ log Q y + + C 


Temperature Distributions with Frictional Heating 

In reference 1 it is shown that the effect of compressibility or 
frictional heating on the generalized velocity distribution is small; it 
was therefore neglected in the preceding section (assumption 6) . In 
this section, the effect of frictional heating on temperature distri- 
butions is investigated and the relation between the velocity and temper- 
ature distributions for a “Prandtl number of 1 is determined. The effect 
of frictional heating in a boundary layer has been previously investi- 
gated by Kalikhman (reference 3). 

In a gas flowing at high velocity, there is, in addition to heat , 
transferred from the wall, heat generated in the fluid by friction. It 
is known that the rate of conversion of flow energy into turbulent and 
heat energy per unit volume is T du/dy. In the present analysis, the 
assumption is made that all the flow energy converted to turbulent energy 
at a point is transformed into heat at the same point, so that the rate 
of generation of heat per unit volume is (t du/dy) /J. For flow over a 
flat surface, this expression is equal to the rate of change of heat 
transfer perpendicular to the surface with distance from the surface, or 
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This equation is applied to flow in a tube inasmuch as for the region 
close to the wall, where the greatest changes of velocity and tempera- 
ture with distance occur, the wall can be considered flat. As in the 
case of velocity distributions, the shear stress is assumed constant 
across the tube (assumption 3) and equation (25) becomes, on integration, 


r u 


pi 


1 

J T o 


du = 


dq 


^10 


or 

<1 “ «o + J T o u (26) 


In obtaining equation (8) in the preceding section, the frictional heat- 
ing was neglected and the heat transfer was considered constant across 
the tube* If the frictional heating is not neglected, equation (8) can 
be written as 


q _ _ dt 

c pS T 0 " 4u 


On substituting the value for q from equation (26) in this equation, 
there results ' 


% + V _ - dt 

c pS T 0 J C " gT n " du 


p° 0 


When this equation is written in dimensionless form, 

a(t/t 0 ) 


3 + 2au + 


du 1 


(27) 


where a is the compressibility parameter, which was used in reference 1. 
The total temperature is defined as 


2gJc p 


T - t + 
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or 



Differentiation gives 




+ 2au + du + 


(28) 


Substituting the value for d(t/t 0 ) from this equation into equation (27) 
gives 



Integration results in 



or . 


0u + = 1 



If T + is defined as 


T + 


(T q - T) c p gT Q 

V T °^ p o 



( 29 ) 


(30) 


equation (29) becomes 


u + = T + 


(31) 
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where T + will "be called the total -temperature parameter. 

Equations (30) and (31) indicate that the velocity and total -temperature 
profiles for compressible flow with heat transfer are similar for a 
Prandtl number of 1. This conclusion was also reached by Kalikhman 
(reference 3). 


Equations for Nusselt Number, Reynolds Number,, 
and Friction Factor 


Bulk temperature . - For obtaining equations for the Nusselt number, 
Reynolds number, and friction factor in terms of the dimensionless para- 
meters used in the preceding sections, an expression for the bulk total, 
temperature is first obtained. The bulk total temperature at a cross 
section in the tube is defined as the temperature the fluid would assume 
in a mixing tank placed immediately downstream of the cross section. 
Thus, the bulk total temperature is given by 


PA 


T v = 


TpudA 


PA 


pudA 


or 


PI Q - 


Tpu (r 0 - y) dy 


nr. 


pu (r Q - y) dy 


(32) 


Equation (32) can be written in dimensionless form as 


t 
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This equation differs from equation (10) because the effects of compressi- 
bility were neglected in obtaining equation (10). Substituting 
equations (31) and (34) into equation (33) gives T^+ for a Eramdtl 
number of 1: 
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Pr 0 + 2 

u +2 (r 0 + - y + ) 


V 


0 


1 - Pu + - au +£ 


dy + 


Pr 0 + 


0 


u + (r Q + - y + ) 

7 ay+ - 

1 - Pu + - au + 


(35) 


Nusselt number. - The heat -transfer coefficient h is defined as 



(36) 


and the Nusselt number with the thermal conductivity based on the wall 
temperature is ' 


Nu 

0 



(37) 


By use of the definitions of T^, r Q + , Pr, and h, the expression for 
the Nusselt number can be written 



2 r 0 + Pr 

V 


(38) 


The Nusselt number with the thermal conductivity based on the wall 
temperature is given by equation (38) . It may be desired to base the 
conductivity on another temperature, perhaps the static bulk temperature 
in order to do so, the ratio of bulk-to-wali static temperature can be 
obtained from 
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rr 0 + 


~ — u+(r. + - y + ) dy + 




t 0 p 0 ° 


Or, 


0 


£- u+ (r n + - y + ) dy + 


Substituting 



and equation (34), this expression becomes 


O r. 


u + (r Q + - y + ) dy + 


b U 


(39) 


or. 


u + (r 0 + - y + ) dy + 


0 1 - 3u + - au +£ 


With the assumptions of constant Prandtl number and constant specific 
heat, the law for the variation with temperature of the conductivity 
must be the same as the law for the variation with temperature of the 
viscosity or 


\ /yf 

k o yo/ 


(40) 


The conductivity at any temperature between the wall and bulk temperature 
can be found in a similar manner. 
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Reynolds number . - The Reynolds number with fluid properties 
evaluated at the vail temperature is defined as 


«e oi V£ 

° p 0 


(41) 


where 


nr r 


*b = 


nr r 


0 


u(r« - y) dy 


(r n - y) ay 


This equation, written in dimensionless form with the denominator 
integrated becomes 


nr o + 


*b 


V T o/ p O 


= V - 


_ + 
r 0 


u + (r Q + - y + ) dy + 


JO 


(42) 


Substituting equation (42) in equation (41) and making use of the 
definition off r Q + give 

»«0 - =W t«) 

The Reynolds number with the fluid properties evaluated at the bulk 
temperature is, from equations (41) and (43), 

<44) 

where %/fci, is calculated from equation (39). 
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Friction factor . - The friction factor with the fluid properties 
evaluated at the wall temperature is defined as 


D(dp/dx) fr _ 2 t q 
2 PQ u h 2 PO u b 2 


(45) 


This equation becomes, in dimensionless form, 



(46) 


where %+ is calculated from equation (42). Another expression for 
f Q can be obtained from equations (38), (43), and (46) if it is assumed 
that u + = T + (Reynolds analogy). This expression is 




Re 


0 


(47) 


From equations (45) and (46), it is evident that the friction factor 
with properties based on the bulk temperature is 



(48) 


Expressions for Nusselt number, Reynolds number, and friction factor 
based on any temperature between the surface and bulk temperature can 
also be obtained by using, the ratio of bulk-to-wall temperature from 
equation (39). 


RESULTS AND DISCUSSION 
Velocity Distributions 

The equations for velocity distributions close to a wall and at a 
distance from a wall for flow with heat transfer (equations (19) and 
(24)) are plotted in figure 1. In this figure, positive values of the 
heat -transfer parameter p indicate heat addition to the gas, whereas 
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negative values indicate heat extraction or cooling. The values for 
the constants (n = 0.109 and k= 0.36), vhich were found in reference 1 
from the experimental data for flow without heat transfer, are used for 
plotting the equations. In reference 1, it was found that the equation 
derived for flow close to a wall without heat transfer agrees closely 
with the data for 0 < y + < 26 and that the equation derived for flow 
at a distance from a wall without heat transfer fits the data for 
y+ > 26. For plotting the present equations for flow with heat .transfer, 
the same limits of applicability for the equations for flow close to a 
wall and at a distance from a wall are used. It can he seen from the 
figure that the exact point of intersection of the curves representing 
the two equations is not critical, especially for high positive values 
of 6, inasmuch as the slopes of the two curves at their intersection do 
not differ greatly. The constant of integration C-^ in equation (24) 

was found for each value of p by substituting at y + = 26 the value 
of u + obtained from equation (19). 

For plotting equation (19), the value of the integral was numeri- 
cally calculated by Simpson's rule. The exponent d in the equation 
was found from viscosity data for air to have an average value of 0.68 
for temperatures between 0° and 2000° F. Although this value was 
specifically obtained for air, the values of d for most other common 
gases do not vary greatly from this value, so that the curves plotted 
in figure 1 should be applicable to most common gases. 

The plots of the equations in figure 1 indicate that increasing 
the rate of heat addition, as is evidenced by increasing the heat- 
transfer parameter p, causes a considerable decrease in values of u + 
for constant y + , as compared with those for adiabatic flow, whereas the 
opposite effect is indicated for heat extraction. A flattening of the 
velocity profile in the central portion of the tube is indicated for 
high rates of heat addition and a peaking of the profile at the center cf 
of the tube is indicated for heat extraction. 

These changes in profile shape are opposite to those for laminar 
flow of a gas ; the explanation is that for turbulent flow, the turbulent 
shear stress, which is governed by the density, is more important in the 
central portion of the tube than is the viscous shear stress. The vari- 
ation with temperature of the density is opposite to that of viscosity, 
so that the changes in profile shape in the. center of the tube for 
turbulent flow caused by heat addition or extraction are opposite to 
those for laminar flow. 
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Nusselt Numbers 

In figures 2 to 5, predicted Nusselt numbers are plotted against 
Reynolds numbers for a fluid with a Prandtl number of 1 for various 
values of the heat-transfer parameter p and the compressibility 
parameter a. In figure 2, the Nusselt numbers and Reynolds numbers 
are plotted with the gas properties (density, viscosity, and thermal 
conductivity) evaluated at the bulk temperature. The Nusselt numbers 
for this case are obtained by dividing equation (38) by the expression 
for kjj/kg given by equation (40); T^+ is calculated from 

equation (35) and t^/tg is calculated from equation (39), where the 

values of the integrals in these equations are numerically calculated 
and r Q +, a, and p are held constant for each integration. The 
relations between u + and y + required in equations (35) and (39) are 
obtained from the plots of equations (19) and (24) in figure 1. The 
Reynolds numbers are obtained from equation (44), where u^+ is calcu- 
lated from equation (42) and t^/to is calculated from equation (39) . 

The exponent d is taken as 0.68 in both equations (40) and (44) because 
the viscosity and the thermal conductivity are assumed to follow the 
same law for variation with temperature. For obtaining various values 
of Nusselt number and Reynolds number for given values of a and p, 
the tube -radius parameter r Q +. is allowed to vary arbitrarily. 

The curves plotted in figure 2, which sire for the case where a 
equals 0 or where compressibility effects caused by high velocities are 
neglected, indicate a decrease in Nusselt number with increasing positive 
values of the heat -transfer parameter p. Increasing values of the heat- 
transfer parameter p correspond, in general, to increasing values of 
the ratio of wall-to-bulk temperature to/t-^; the curves shown correspond 
to a range of ratios of wall-to-bulk temperature of from 1 to about 7 for 
heating and from 1 to about 0.4 for cooling.' The same trends for heating 
have been experimentally observed to determine over -all heat -transfer 
coefficients in tubes with a range of ratios of surface -to-bulk tempera- 
ture of from 1 to about 3^- (references 4 to 8). The dashed line in 

figure 2 represents the mean of experimental data for moderate tempera- 
ture differences (reference 9) and is seen to agree, in general, with 
the predicted line for p = 0. 

For heat addition to the gas, figure 3(a) shows predicted Nusselt 
numbers plotted against Reynolds numbers for two values of the heat- 
transfer parameter p and for a = 0 with the gas properties in the 
Nusselt and Reynolds numbers evaluated at a film temperature between the 
wall and bulk temperatures defined by 
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*0.4 ^ °* 4 (t 0 ' *b) + *b 


(49) 


or 



0.4 


+ 



For obtaining the values for Nusselt and Reynolds numbers for figure 3(a), 
the same equations were used as were used for figure 2 except that t^/tQ 
was replaced by tQ^/tQ in the equations. The plot in figure 3(a) 
indicates that the effects of p, or of ratio of wall-to-bulk temperature, 
is practically eliminated for heat addition when the gas properties are 
evaluated at the film temperature as defined. In the experimental work 
described in references 4 to 8, it was found that the effects of ratio of 
wall-to-bulk temperature on the Nusselt numbers were, practically eliminated 
when the fluid properties were evaluated at either the wall or the average 
film temperature. In the present analysis, however, where the m a ximum 
ratio of wall-to-bulk temperature was about twice that in the experimental 
work, separation of the curves was found when the properties were evaluated 
at the wall temperature . 

The curve from figure 3(a), together with- experimental data from 
reference 8, is shown in figure 3(b)'. The data points were calculated in 
reference 8 by evaluating the gas properties at the average of the wall 
and bulk temperatures; substantially the same results are obtained if the 
properties are evaluated at t Q>4 for the range of ratios of wall-to- 
bulk temperature used in the tests. The data, which were obtained for 
air, were corrected to a Prandtl number of 1 by dividing the Nusselt 
numbers by the empirical correction Pr0*4. The data are seen to agree 
substantially with the predicted line for heating for Reynolds numbers 
greater than 10,000. The deviation of the data from the predicted line at 
low Reynolds numbers is probably caused by transitional effects in the fluid 
flow that were not considered in the analysis. It should be mentioned 
that the analytical and experimental results are not strictly comparable 
because the analytical results are for local heat -transfer coefficients 
at a cross section in the tube, whereas the experimental results give 
average heat-transfer coefficients for a whole tube. 

% 

Variations of predicted Nusselt numbers with Reynolds numbers for 
flow with heat extraction or cooling are shown in figure 3(c) with the 
properties evaluated at a film temperature defined by 


*0.6 - 0,6 (*0 " *b) + *b 
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It Is seen that use of this temperature eliminates the effects of 0 
on the Nusselt numbers for heat extraction. 

A correlation for both heating and cooling in figure 3(d) shows 
Nusselt numbers with the gas properties evaluated at the average of 
the wall and bulk temperatures to. 5 (This temperature is also the 

average of tQ ^ and tg g.). By use of this temperature for evaluat- 
ing the fluid properties, the maximum separation of the curves la found 
to be about 9 percent. 

In figures 4 and 5 are shown predicted Nusselt numbers plotted 
against Reynolds numbers for flow with heat addition where the compressi- 
bility effects associated with hlgi subsonic velocities are considerable 
(a = 0.00025). Far plotting these curves, the velocity-distribution 
equations (equations (19) and (24)), which were derived far a s 0, were 
used because the effects of compressibility on velocity distributions 
were found in reference 1 to be small for subsonic flow. As In the case 
of flow for a = 0, the use of the static film temperature i, defined 
by equation (49), for evaluating the gas properties practically eliminates 
the effect of the heat-transfer parameter 0 on the Nusselt numbers for 
heat addition. Part of the slight difference between the plots in fig- 
ures 3 and 5 is probably caused by assuming that a = 0 for the velocity- 
distribution equations for. both plots. 

The heat -transfer coefficients used in figures 2 to 5 are defined 
by equation (36) and are based on the difference between the wall and > 
bulk total temperatures. Although strictly speaking equation. (36) should 
be used only for fluids with a Prandtl number of 1, it can be used for 
real gases except for cases where the temperature differences are very 
small and the Mach numbers are high (reference 10) . 


Friction Factors 

Predicted friction factors plotted against Reynolds numbers for 
various values of the heat -transfer parameter 0, with a = 0, are shown 
in figures 6 and 7. In figure 6, the friction factors are calculated 
from equation (48) where u^+ and t^/tQ are calculated as before from 

equations (42) and" (39). As in the case of Nusselt numbers, a decrease 
in friction factor is indicated, for- high positive values of 0 or of 
ratio of wall-to-bulk temperature when the gas properties in the friction 
factors and Reynolds numbers are evaluated at the bulk temperature. 

These trends have also been experimentally observed (reference 5). The 
dashed line in figure 6 represents the mean of experimental data for 
isothermal flow (reference 11) and agrees in general with the pr edic ted 
line for Reynolds numbers above 10,000. 
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Friction factors- are plotted against Reynolds numbers in figure 7(a) 
for heat addition with the gas properties evaluated at the film tempera- 
ture defined by equation (49), indicating that, as in the case of 
Nus se It numbers, evaluating the gas properties at that temperature 
practically eliminates the effects of p or of ratio of wall -to -bulk • 
temperature for flow with heat addition. In figure 7(b), which is 
plotted for flow with heat extraction, it is indicated that the use of 
^0.6 for evaluating the gas properties eliminates the effects of (3. 
or of variable gas properties on the friction factors. for heat extraction. 
The friction-factor correlation for both heating and cooling of fig- 
ure 7(c) indicates that, as in the case of Nusselt numbers, the 
separation of the curves for various values of p is about 9 percent 
when the gas properties are evaluated at the average of the wall and 
bulk temperatures. 


SUMMARY OF RESULTS 

From an analytical investigation of fluids flowing through smooth 
tubes with heat transfer, the following results were obtained; 

1. The equations for the prediction of velocity distributions for 
fully developed turbulent flow of gases in smooth tubes with heat transfer 
indicated that variation of fluid properties across the tube causes a 
flattening of the velocity profile for heat addition and a peaking of the 
profile at the center of the tube for heat extraction. 

2. The investigation of the effect of frictional heating on tempera- 
ture distributions indicated that, for a fluid with a Prandtl number of 1, 
the velocity and total -temperature profiles are similar. 

3. The equations for Nusselt number, Reynolds number, and friction 
factor for a fluid with a Prandtl number of 1 for the case where the 
variation of fluid properties across the tube is large predicted trends 
that agree with those experimentally obtained for Reynolds numbers ' 
greater than 10,000. 

4. The analysis indicated that the effect of variation of fluid 
properties across the tube on the Nusselt number and the friction-factor 
correlations for both heating and cooling can be eliminated by evaluating 
the fluid properties at a temperature close to the average of the wall 
and bulk temperatures . 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, July 28, 1950. 
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APPENDIX - SYMBOLS 

The following symbols are vised in the report: 

. ' cross-sectional area based on inside diameter of 

tube, sq ft 


C, Cj 

C P 

D 

d 

g 

h 


constants of integration 

specific heat of fluid at constant pressure, 
Btu/(lb)(°F) 

inside diameter of tube, ft 1 

exponent that describes variation of viscosity of 
fluid with temperature 

acceleration due to gravity, 32.2 ft/sec 2 

<10 


heat-transfer coefficient 

<1q 


T 0 - 


for compressible 


fluid or 


*0 " t b 


for incompressible fluid, 


Btu/(sec)(sq ft)(°F) 

J mechanical equivalent of heat, 778 ft-lb/Btu 

k thermal conductivity of fluid,. (Btu)(ft)/ 

(sec)(sq ft)(°F) 

k, thermal conductivity of fluid evaluated at t. , 

’ (Btu)(ft)/(sec) (sq ft)(°F) 

kg thermal conductivity of fluid evaluated at wall, 

(Btu)(ft)/(sec) (sq ft)(°F) 

k 0.4> k 0.5> k 0.6 thermal conductivity of fluid evaluated at tg, t Q 4 , 

tg >5 , and tg >6 , respectively, (Btu)(ft)/ 

(sec)(sq ft)(°F) 


n 


constant 

total rate of heat transfer toward tube center per 
unit area, Btu/(sec)(sq ft) 
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<*t 


In . 


T 0 or * 0 ' 


0.4 


'0.5 


' 0.6 


u 


u b 

y 

€ 

€ h 

n 

^b 

^0 


rate of turbulent heat transfer toward tube center per unit 

area, Btu/(sec)(sq ft) 

\ 

rate of viscous heat transfer toward tube center per unit 
, area, Btu/(sec)(sq ft) 

total rate of heat transfer at wall toward tube center per 
unit area, Btu/(sec)(sq ft) 

inside tube radius, ft 

total temperature, °R 

bult or average total temperature of fluid at cross section . 
of tube,- °R 

\ 

absolute wall temperature, °B 
absolute static temperature, °R 

bulk: or average static temperature of fluid at cross section 
of tube, °B 

film temperature, 0.4 (t Q - t-j,) + t^, °R 

film temperature, 0.5 (tp - t^) + tjj, °R 

film temperature, 0.6 (t^ - t^) + t^, °R 

time average velocity parallel to axis of tube, ft/sec 

bulk or average velocity at cross section of tube, ft/sec 

distance from tube wall, ft 

coefficient of eddy diffusivity for momentum, sq ft/sec 
coefficient of eddy diffusivity for heat, sq ft/sec 
Karmen constant 

absolute viscosity of fluid, (lb)(sec)/(sq ft) 

absolute viscosity of fluid evaluated at t^, (lb)(sec)/(sq ft) 

absolute viscosity of fluid at wall, (lb)(sec)/(sq ft) 
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Uq.4> Pq. 5 > ^0.6 afcs °l u ' t '© viscosity of fluid evaluated at t 0> 4 , 

*0.5* 311(1 t 0.6> respectively, (lb) (sec)/(sq ft) 

p mass density, lb-sec^/ft^ 

P-u bulk or average density at cross section of .tube, 

lb-sec 2 /ft^ 

Pq mass density of fluid at wall, lb-sec 2 /ft 4 

Pq , 4 * Pq. 5 ^ Po. 6 densit y of fluid evaluated at t 0>4 , t Q>5 , and 

to > 6; respectively, lb-sec^/ft^ 

I 

T shear stress in fluid, lb/sq. ft 

T 0 shear stress in fluid at wall, lb/sq ft 

Subscript : 

fr on friction pressure gradient . 

Dimensionless groups: 

a compressibility or frictional heating parameter, 

T 0 

2 6 Jc p T 0 p 0 

(3 heat -transfer parameter, _ , 

°p 8 Vo 

D(dp/ta). 2 t. 

f_ friction factor, — = 

2o oV p oV 

f^ friction factor with density evaluated at t-^, 

D(dp/da) fr 2 t q 

2 Pb u b 2 Pb u b 2 

t- m9 r-, f rt ~ friction factor with density evaluated at t n t n 
0.4- 7 0.5 7 0.6 U.4 7 O.o 7 

and tQ g, respectively 
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Nu b 

< Nusselt number with thermal conductivity evaluated 

a + f hD 

az 7T- 

*b 

Nii 0 

Nusselt number with thermal conductivity evaluated ' 
at tg 

Wu 0.4> Nu 0.5> Nu 0.6 

Nusselt number with thermal conductivity evaluated 


*o. 4 , to. 5 , and to. 6 , respectively , 

Pr 

c #g 

Prandtl number, — 

k 

Re b 

Reynolds number with density and viscosity evaluated 
at t 'W > 

a * V p b 

Ee 0 

Reynolds number with density and viscosity evaluated 
at tg 

t 

Re 0 . 4 > Re 0 #5 , Re 0>6 

Reynolds number with density and viscosity evaluated 
at *0.4' *0.5' a " 1 * 0 . 6 ' res P ectivel y 

r + 
0 

tube -radius parameter, 9^2. r n 

^o/ p o 

T + 

(T 0 - T)c_gr 0 

total -temperature parameter. — t — — 

^0 W p 0 

m + 

a d 

bulk total -temperature parameter, i ( 1 - _Jl\ 

P V V 

U + - 

velocity parameter. u 

yV p o 

V. 

bulk velocity parameter, Ulj 

l/V p o 

y + 

V T o^ p o 

vall-distance parameter, 7 v 

P 0' P 0 


y 
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q H/CU ■= q nn 


Figure 2 . - Variation of Nusselt number with Reynolds number for flow of gases with gas properties evaluated at static bulk temperature. Compressibility 
parameter a, 0; Prandtl number, 1. Positive and negative values of 3 Indicate heat addition and extraction, respectively. 







Figure 3. - Continued, Variation of Nusaelt number with Reynolds number for flow of gases. Compressibility 
, parameter o, 0; Prandtl number, 1. 




Figure 3. - Continued. Variation of Nusaelt number with Reynolds number for flow of gases. Compressibility parameter a, 0; Prandtl 

♦ number, 1 . 
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(d) Evaluation of gas properties at to . 5 3 0.5(to~tb)+ t^; heat addition and extraction. 




Figure 4* - Variation of Nusselt number with Keynolds number for flow of gases at high velocities with heat addition with gas properties evaluated at static 

bulk temperature. Compressibility parameter a, 0.00025; Prandtl number, 1, 
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Figure 5. - Variation of Nusaelt number with Reynolds number for flow of gases at high velocities with heat addition with gas properties evaluated 

at to , 4 “ 0.4(tQ-tb) ♦ Compressibility parameter a, 0.00025; Prandtl number, 1 . 





Figure 6. - Variation of friction factor with Reynolds number for flow of gases with gas properties evaluated at static bulk temperature. Compressibility 
parameter a, 0; Prandtl number, 1. Positive and negative values of g Indicate heat addition and extraction, respectively. 










Figure 7. - Continued. Variation of friction factor with Reynolds number for flow of gases. Compressibility parameter 

Prandtl number, 1 . 





Figure 7. - Concluded- Variation of friction factor with Reynolds number for flow of gases. Compressibility parameter c # 0; Prandtl number. 



